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ABSTRACT 

We present the occultation observation of compact radio source B0019— 000 through the plasma tail of comet Schwassmann- 
Wachmann 3-B. The observation was made with the Ooty Radio Telescope at 326.5 MHz on 2006 May 26, when the plasma tail 
of the comet was in front of this source. The scintillation was found to be increased significantly for the target source compared to 
that of a control source. The intensity fluctuation power spectra show both steepening at high spatial scales and excess power at 
low spatial scales. This observation can be attributed to the turbulence in the comet plasma tail. A two-regime plasma turbulence 
can explain the time-evolution of the power spectrum during the occultation observation. 

Subject headings: comets: general — comets: individual (Schwassmann-Wachmann 3) — interplanetary 
medium — occultations — turbulence 
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1. INTRODUCTION 

The interplanetary scintillation (IPS) technique is a useful, 
convenient, remote-sensing method to probe th e solar wind 
dHewish. Scott. & WillslfT964l:lManoharanlll993h . It exploits 
the scattering of radio waves from a compact radio source by 
the density irregularities in the solar wind and provides infor- 
mation on the speed and turbulence spectrum of the solar wind 
plasma at the region of cl osest solar offset o f the line of sight 
to the radio source (e.g., Manoharanl l 1993b . When a comet 
plasma tail intersects the radio path, the electron density ir- 
regularities associated with the comet can cause an enhance- 
ment in the level of scintillation. Several attempts have been 
made to estimate the turbulence associated with the cometary 
plasma by observing the occultation of radio source by the tail 
of co met (lA nanthakrishn an et aD 1 1 975L 119871: lAlurkur et alJ 
19861 U989hlSlee et alJll987t lHaiivassil iou. & Duffett-Smithl 
1987b . However, in order to confirm the effect of the cometary 



plasma as well as to assess the contribution to the level of tur- 
bulence imposed by any solar wind transients or disturbances, 
it is essential to simultaneously monitor a region just outside 
th e tail of the comet during such observations. 

(lAna nthakrish nan et alJ 1 19871) reported a significant in- 
crease in the plasma turbulence as the tail of comet Halley ap- 
proached a radio source. But, a rather similar enhancement in 
scintillation was also observed towards c ontrol sources, which 
were located outside the comet tail ( lAnanth akrishn an" et alJ 
1 19871) . This emphasized the need for near-simultaneous mon- 
itoring of IPS through the tail of the comet as well as a region 
outside the comet tail. This study attributed an insignificant 
increase of turb ulence associated with the cometary plasma. 
In a later work, ISlee et all d 1990b adopted the idea of near- 
simultaneous observation of control source during cometary 
occultation observation and reported a significant increase of 
the scintillation of the target source compared to that of the 
control source. In this Letter we report the occultation of the 
compact radio source B0019— 000 by the plasma tail of comet 
Schwassmann-Wachmann 3-B. The results are substantiate by 
the observations of the IPS of a control source outside the tail. 
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Systematic changes observed in the IPS spectra as the source 
approached the comet tail suggest a steeper turbulence spec- 
trum for the cometary plasma at scales >500 km and exces- 
sive turbulence at the small-scale part (i.e., <50 km) of the 
spectrum. Below, in §2, we give the details of our observa- 
tions and the results. We compare the results with previously 
reported results and discuss the actual astrophysical implica- 
tion of the results in §3. Finally, we summarize and present 
our conclusions in §4. 

2. OBSERVATIONS AND RESULTS 

The comet Schwassmann-Wachmann 3 (also known as 
73P/Schwassmann-Wachmann) is a periodic comet. This 
Jupiter-family comet passed within ^0.1 AU of the 
Earth around the middle of 2006. It is to be noted 
that the comet was going through breaking and disinte- 
gration dBoehnhardt. & Kaufll 119951: ICrovisier et all [1996; 
Boehnhardt et a 1. 2002), which would likely to increase the 
level of turbulence. The occultation observation reported 
here was therefore aimed at estimating the turbulence spec- 
trum of density irregularities and associated scale sizes of the 
cometary plasma. 

Compact sources lying close to the path of the comet were 
chosen from the Molonglo Catalogue (Large et al. 1983). On 
2006 May 26, one of the fragments (labeled 73P-B, the sec- 
ond brightest fragment at the time of observation) aligned 
along the line of sight to the radio source BOO 19— 000, which 
has a compact component (~ 50 mas) containing > 75% of 
the total flux density, S 40 8=2.99 ±0.14 Jy. The line of sight 
along the radio source B0019-073 (S 408 = 1.83±0.07 Jy), 
was monitored as the control region. The systematic moni- 
toring of these sources on the day of occultation as well as 
for the next two consecutive days provide the spectrum of the 
density irregularities and dominant scale sizes of the cometary 
plasma. This Target of Opportunity observation was carried 
out using the Ooty Radio Telescope (ORT) at the Radio As- 
tronomy Centre, Ooty dSwarup et alj|1971h . The ORT, an 
equatorially mounted parabolic cylindrical antenna, is 530 m 
long in the north-south direction and 30 m wide in east-west 
direction. It operates at 326.5 MHz and can effectively probe 
solar wind in the heliocentric distance range of 0.05 -1.0 AU. 

On 2006 May 26 - 28, both the target source (B0019-000) 
and the control source (B0019— 073) were observed. The co- 
ordinates and flux densities of these sources are given in Ta- 
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TABLE 1 
Radio sources parameters 



IAU Name 


Ecliptic coordinates 


^408 


B1950 J2000 


J2000 (deg) 


(Jy) 


0019-000 0022+002 


5.245 - 1.998 


2.99±0.14 


0019-073 0022-070 


2.320 - 8.685 


1.83±0.07 



a From Molonglo Reference Catalogue (Large et al. 1981) 



TABLE 2 

Comet 73P/Schwassmann-Wachmann3-B parameters' 1 



Date 


Nuclear ecliptic 


Heliocentric 


Topocentric 


Elongation 


at 


coordinates 


Distance 


Distance 


Leading 


0800 UT 


(dcg) (deg) 


(AU) 


(AU) 


(dcg) 


May 26 


251.15 -0.25 


0.957 


0.1216 


59.62 


May 27 


252.48 -0.52 


0.955 


0.1287 


59.60 


May 28 


253.82 -0.79 


0.952 


0.1360 


59.70 



a From JPL solar system dynamics group 



ble ([T). The orbital parameters of the comet Schwassmann- 
Wachmann 3-B are summarized in Table ©. Fig. (Q]) shows 
the position of the nucleus of the comet on 2006 May 26 - 
28 with respect to the position of both the target and con- 
trol sources, projected on the plane of the sky. Fig. (ffji dis- 
plays the comet trajectory on 2006 May 26 over-plotted on 
the 1.4 GHz image of the field centered at the target source 
B0019-000 taken f rom the NRAO VLA Sky Survey (NVSS; 
ICondon et al.lll998h . The positions of the comet on these days 
are taken from the online high-precision ephemerides pro- 
vided by Jet Propulsion Laboratory (JPL) solar system dy- 
namics group. The target source was observed during occul- 
tation on 2006 May 26 from 07 1 3 to 0804 UT (when the target 
hour angle reached the west limit of the ORT) at 326.5 MHz 
with a total bandwidth of 4.0 MHz and a time resolution of 20 
ms. During this span, the control source was observed for a 
short period to understand the solar wind characteristics and 
to distinguish between any changes caused by the turbulence 
in the tail of the comet and the solar wind. Note that the con- 
trol source is sufficiently close to the target source. So, these 
two lines of sight are expected to have similar solar wind char- 
acteristics. A nearby cold region of the sky was also observed 
to get the "off-source" reference. Both the target and the con- 
trol sources were observed again on 2006 May 27 and 28. 

A detailed account of the IPS observ ation using the ORT 
and data analysis procedure is given in iManoharanl (120031) . 
In brief, the Fourier transformation and spectrum computa- 
tion from the recorded intensity scintillation measurements 
are performed to get the power spectrum of the intensity fluc- 
tuations to bring out the statistical properties of the interven- 
ing medium. A 2048-point fast Fourier transformation yields 
a frequency resolution of «0.025 Hz. The spectrum is cor- 
rected for the off-source noise and the averaging of the ad- 
jacent spectral points gives a uniform noise up to the cutoff 
frequency (where the power drops to noise level). The area 
under the power spectrum is the measure of the scintillation 
index. Since the present measurements are made in the weak 
scattering region, the shape of the spectrum and the scintilla- 
tion index are directly related to the turbulence spectrum of 
the interve ning plasma and th e level of plasma turbulence, re- 
spectively ( Ma noharanl 1993b . 

As illustrated in Fig. ©, significant changes are evident 
in the shape of the power spectrum as the comet nucleus ap- 
proached the line of sight to the radio source. The time of ob- 
servation (in UT) for each spectrum is given on the top right of 
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Fig. 1 . — Path of the comet with respect to the radio sources (both 
the target and the control source). The position of the comet nucleus 
on2006 May 26 - 28 is plotted. Part of the trajectory on 2006 May 
26 is shown in the inset, where the positions close to the time of 
observation are marked. 



NVSS: 0019-000 (levs = + /- 1 , 1 .4,2. 2.8,4. ..mjy/b) 




Right Ascension(J2000.0) 

Fig. 2. — NVSS L-band image of the field centered at the target 
source 0019—000. Part of the trajectory on 2006 May 26 is over- 
plotted with the position at 0600, 0700, 0800, and 0900 UT (top 
right to bottom left) marked with crosses. 

each panel. The details of the implication of the power spec- 
trum shape and the interpretati on of scale size from IPS spec- 
tral va riati on are described in |Manoharan. Koiima, & Misawa 
d 19941) and lManoharan et al.1 (1200(1 7 Panels a - c in Fig. © 
show the change of the power spectra during occultation on 
2006 May 26. Panels d and e are the power spectra from the 
follow-up observations on 2006 May 27 and 28, respectively. 
The power spectra for the control source taken on May 26 
(near simultaneous to occultation observation), 27, and 28 are 
shown in panels f — h, respectively. The spectrum towards the 
source obtained from the first half of the observation (0713 - 
0740 UT) on May 26 has no significant change. But spec- 
tra from the later part of the observation (0743 - 0804 UT) 
shows steepening at lower spatial frequency. The scale size 
a = v/(2%f), where / is the temporal frequency and v is 
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Fig. 3. — Power spectra of the intensity fluctuation. The horizontal axis is the spectral frequency (Hz) and the vertical axis is power (dB). 
Date, source name, and the time of observation (UT) are given at the top of the each panels. The dashed lines indicate the noise power level. 
Note the steepening of the spectrum at high spatial scales (panel b and c) and significant excess power at low spatial scales (panel c) compared 
to the rest of the spectra from follow-up and control observations (See Sj2]for details). 



0.9 



X 


~o 
_c 

c 

o 



o 

00 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 
0019-000 




TABLE 3 

Parameters for the plasma tail on May 26 



Time 


Flux a 


Scintillation 


AS b 


AS 2 


(*no/n 2 ) c 


(UT) 


(Jy) 


index 


(mJy) 


(AS 2 ) 


(*nc/n 2 ) s 


0725 


3.1 


0.068 


210.8 


1.327 


17.8 


0750 


3.0 


0.071 


213.0 


1.355 


19.3 


0800 


2.7 


0.089 


240.3 


1.724 


39.4 



250 



a Flux density measurement is limited by the ORT confusion. 
b AS is not limited by confusion and is much accurate. 

L, the scintillating flux density variance AS 2 oc $ nc L. From 
this, we can relate the change of measured scintillation index 
to the physical parameters of solar wind and cometary plasma, 
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Fig. 4. — Scintillation index enhancement in occultation. The solid 
line represents the long-time average of the scintillation index as a 
function of heliocentric distance of the point of closest approach of 
the line of sight towards 0019—000 obtained from scintillation moni- 
toring for many years from the ORT. The triangles are measurements 
during occultation. The part of the plot of our interest is magnified 
in inset. 

the solar wind speed. Considering a typical solar wind speed 
of ^400 km s _1 , this translates to spatial scales of > 500 
km. Excess power is also seen at low spatial scales (< 50 
km) on the spectrum obtained from the end part of the obser- 
vation (0754 - 0804 UT). The shape of the power spectra of 
the control observations shows normal quiet solar wind con- 
dition. The changes in the spectra during the occultation are 
attributed to the turbulence in the comet plasma tail. 

We find that the scintillation index of the control source was 
normal during these observations and also during the control 
observations on the following days, implying normal quiet so- 
lar wind conditions. But we have noticed significant changes 
in the scintillation index for the target source. The scintillation 
index increased, in about 1 hr time during the course of our 
observation, from 0.068 to 0.089. It is clear from the power 
spectra that the steepening at the high spatial scales is the 
dominant source of this enhancement. Here we consider the 
ratio of the level of the electron density spectrum $ nc to the 
mean squared electron density (n 2 ) to be a measure of plasma 
turbulence. For an intervening medium of effective thickness 



AS 2 {<P nc /n 2 ) s n 2 cs L s + (<t> nc /n 2 c ) c n 2 cc L c 
{AS 2 ) ($ ae /nl) s nlL s 



(1) 



where the subscripts Y and 'c' refer to solar and cometary, 
respectively, and (AS 2 ) is the estimated variance due to the 
solar wind on that day. 

We have used the data obtained from scintillation moni- 
toring observations from the ORT for many years towards 
0019-000 to get the scintillation index as a function of helio- 
centric dist ance of the point of clo sest approach of the line of 
sight (see Manoharanl ll993L 120 06, for details of scintillation 
index variation with heliocentric distance). During our ob- 
servation the average heliocentric distance was ~ 185.5R©. 
Fig. (01 shows the long-time average of the scintillation index, 
and from this, for normal solar wind conditions, the expected 
scintillation index is 0.061 at 185. 5R© for a flux density of 
3.0 Jy. This corresponds to (AS" 2 ) - 5 = 183 mJy. Please note 
that the peak scintillation index is ~ 0.8 for this source. This 
implies that the source is very compact (< 50 mas) and is an 
ideal source for scintillation measurements. From our occul- 
tation data we can find, using Eqn. (Q~|), the normalized level 
of turbulence ($ nc /?i 2 ) in the comet plasma tail compared to 
that in the solar wind. We assume typical values for L s and L c 
to be 10 s and 10 5 km, respectively. Considering that the radio 
source line of sight probes the plasma tail at about 1.4 x 10 4 
km away from the nucleus and that the angular size of the 
coma was ~ 3', the expected value of L c is ~ 5 x 10 4 km. 
L c = 10 5 km may be taken as an upper limit and, for smaller 
values of L c , the corresponding value of (Q nc /n 2 ) c will in- 
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crease. Average electron number densit y in the cometar y 
plasma is assu med to be n ec = 30 cm " 3 dGrard et al.|[l986h. 
Earlie r reports dAlurkur et al.l 19861; ?; [A nanthakrishnan et al. 
119871) used, following iReme et all ( fl986l) . n ag = 5 cm" 3 for 
normal solar wind condition. We have instead taken the av- 
erage value of number density to be 7 cm' 3 computed from 
the 1 minute OMNI data (from the CDAWeb service of the 
NASA/GSFC Space Physics Data Facility) for the period of 
2006 May 21-31. We summarize the result of plasma turbu- 
lence enhancement in Table ©. 

3. DISCUSSION 

We present the detection of a rare event of the occulta- 
tion of a compact source (B0019— 000) by the plasma tail of 
comet Schwassmann-Wachmann 3-B using scintillation ob- 
servation and the power spectral analysis. It is interesting to 
compare the present results with the earlier reports of both 
detections and non-detections of enhanced scintillation. The 
range of values of the normalized plasma turbulence from 
similar observations spread o ver almost an order of mag- 
nitude (lAlurkuretal.iri986l ?: lAnanthakrishnan et ai]|1987t; 
ISlee et alJI 1990b for a given set of parameters (e.g., L c , n ec ). 
In addition to that, the turbulence level is likely to change sys- 
tematically with the distance from the nucleus of the comet 
along the plasma tail axis. In other words, occultation ob- 
servations probe the plasma at different distances from the 
comet nucleus and hence may reveal a variety of turbulence 
levels. Hence a straightforward quantitative comparison of 
plasma turbulence is not possible. But from earlier results, 
the clear trend emerges that most of the negative results can 
be attributed to a chance of actually probing a lower level of 
turbulence in the occultation events dSlee et alJI 1990b for var- 
ious reasons. The rest of the cases with favorable observing 
geometry have resulted in a positive detection. In the present 
study, given the possibility of enhanced level of turbulence 
because of the ongoing disintegration and the small offset of 
the line of sight from the nucleus of the comet, our result is in 
good agreement with this scenario. 

The change of the intensity fluctuation power spectra has 
contributions in two distinct scales. The steepening of the 
spectrum at low frequency indicates scale sizes in the inter- 
vening plasma at high spatial scales and remarkable turbu- 
lence enhancement associated in this spectral region. For a 
typical solar wind speed of 400 km s _1 , it translates to linear 
scales of > 500 km. At the later part of the observation, at the 
close encounter of the comet nucleus and the line of sight to 
the radio source, we notice the additional significant contribu- 



tion in the power spectrum at the low-frequency region, which 
corresponds to linear scales of < 50 km. It is to be noted that 
the IPS measurements normally show much less turbulence at 
these small scales. The evolution of the spectrum with time 
during the course of the occultation indicates that the plasma 
turbulence level is dominated by small-scale density irregu- 
larities closer to the tail axis. However, it is the large-scale 
irregularities, developed most likely due to the interaction of 
plasma with the ambient solar wind and the evaporation pro- 
cess, that contribute significantly at the outer part of the tail. 
This shows the presence of two main regimes of turbulence in 
the comet plasma tail. This is found to be in ag reement with 
the earlier reported results of ^Slee et ail (1 19901) for the cases 
of comet Halley and comet Wilson, where similar two-regime 
plasma turbulences were detected. 

4. CONCLUSIONS 

We have present the results from the occultation observa- 
tion of a compact radio source B0019— 000 at 326.5 MHz 
through the plasma tail of comet Schwassmann-Wachmann 
3-B. We detect significant enhancement in the scintillation in- 
dex and change in the intensity fluctuation power spectrum 
during the occultation period. The systematic study of the 
evolution of power spectrum clearly reveals the dominant ef- 
fect of large-scale (> 500 km) density fluctuations and the 
presence of a smaller scale (< 50 km) closer to the tail axis. 
Since the shape of the power spectra of all the control observa- 
tions show normal quiet solar wind condition, these changes 
at the time of occultation are attributed to the plasma tail of 
the comet. This clearly shows the presence of a two-regime 
plasma turbulence with two distinct scales of density irregu- 
larities that may arise from the interaction of solar wind with 
the comet plasma tail. 
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